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■ 1 Introduction 

cr 

The cosmos of general relativity is a pseudo-Riemannian manifold (V, g) of Lorentzian signature 
(—,+,... ,+). The Einstein equations link its Ricci tensor with a phenomenological tensor 
^ , which describes the stresses and energy of the sources. They read: 

Ricci(#) = p, 

that is in local coordinates x x , A = 0, 1, 2, . . . , n, where g = gx^dx^dx 11 (classical physics n = 3), 

R <*P = ~Qx^ Ta P ~ dx^ T > 3X + r «/3 r V ~ r ^ r /3A = Paf3> ( L1 ) 

where the T's are the Christoffel symbols: 

and p is a symmetric 2-tensor given in terms of the stress energy tensor T by, 

P a p = T a/3 ~ ^g a /3 tr T, with tr T = g XfJ, Tx^. 
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Due to the Bianchi identities the left hand side of the Einstein equations satisfies the 
identities, with V a the covariant derivative in the metric g, 

V a (R?e - \g aP R) = 0, R = g^R x » 

The stress energy tensor of the sources satisfies the conservation laws which make the equations 
compatible, 

V a T a/3 = 0. 

In vacuum the stress energy tensor is identically zero. The presence of sources brings up new 
problems specific to various types of sources. 

Since very little is known about the global properties of the universe, it is legitimate to 
study arbitrary manifolds and metrics. Also, in modern attempts to unify all the fundamental 
interactions, manifolds of dimension N greater than four, endowed with metrics with (N — 4)- 
dimensional isometry groups are considered. We will therefore not restrict our study to four- 
dimensional manifolds, when possible. 

Since we will treat only non-quantum fields, it seems that a first problem to look at is 
the problem of classical dynamics, i.e., the problem of evolution of initial data. The Einstein 
equations are a geometric system, invariant by diffeomorphisms of V, the associated isometries 
of g, and transformation of the sources. From the analyst's point of view they constitute, for 
the metric, a system of second order quasi-linear partial differential equations. The system is 
determined because the characteristic determinant is identically zero (a property linked with 
diffeomorphism invariance). It is overdetermined because the Cauchy data are not arbitrary 
(which must be interpreted geometrically). 

To study the Cauchy problem we must split space and time. The evolution will be formu- 
lated for time dependent space tensors. 

2 Moving frame formulas 
2.1 Frame and coframe 

A moving frame in a subset U of a differentiable (n + l)-dimensional manifold V is a set of 
(n+ 1) vector fields on U linearly independent in the tangent space T X V at each point x eU. A 
coframe on U is a set of (n+ 1) 1-forms 9 a linearly independent at each x € U in the dual T*V. 
In the domain U of a chart a coframe is defined by (n + 1) linearly independent differential 
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1- forms, 

9 a = a%dx fi , (2.1) 

with dp functions on U. 

The metric is Lorentzian if the quadratic form is of Lorentzian signature. 

Remark 2.1 The splitting V = A4 xIR, with M an orientable 3-manifold, implies the existence 
of a global coframe (but not of global coordinates!). 

The coframe defined by the 1-forms 6 a is called the natural frame if 9 a = dx a . 

In a general frame the differentials of the 1-forms 9 a do not vanish; they are given by the 

2- forms, 



dd a = C| 7 ^ A # 7 . (2.2) 



The functions (7g on hi are called the structure coefficients of the frame. 
The Pfaff derivative d a of a function on hi is such that, 

df = ^dx a = d a f6 a . (2.3) 

We denote by A with elements Aa of the matrix inverse of a with elements af. It holds that: 



Pfaff derivatives do not commute. One deduces from (|2.2[) and the identity d f = that 



d 2 f = \{d^f - d^dpf - C^d a f}6? A9^ = 0, (2.5) 

hence, 

(d a dp - d p d a )f = ClpOyf. (2.6) 

2.2 Metric 

A metric on hi is a nondegenerate quadratic form of the 6 a, s: 

g = g aP e a e^. (2.7) 

A frame is called orthonormal for the metric g if g a p = ±1. In the case of a Lorentzian metric 
we will denote by 6° the timelike (co)axis and l the space (co)axis, then in an orthonormal 
frame, goo = — 1 and gij = 5ij, the Euclidean metric. 



2.3 Connection 

A linear connection on V permits the definition of an intrinsic derivation of vectors and tensors. 
It is denned in the domain IA by a matrix- valued 1-form u i.e., by a set of matrices u>^ linked 
to functions w« 7 by the identities, 

u% = u^9 a . (2.8) 
The covariant derivative of a vector v with components v a is, 

V a ^ = 9a^+W^«T. (2.9) 

An analogous formula holds for a covariant vector now with a minus sign in front of w^y. 
Definition 2.1 The connection uj is called the Riemannian connection of g if 

• It has vanishing torsion, i.e., 

den + ul p e a AO? = o, (2.io) 

that i^, 

w /3 7 ~~ W 7/3 — u /3 7 - 

• The covariant derivative of the metric is zero, i.e., 

da3/3 7 - ^ay9/3X ~ ^a^X-y = 0- (2.11) 

These two conditions imply by straightforward computation that, 

^ 7 = r£ 7 + /^ a7 , M , (2.12) 

with 

■^(dtixCa-y ~ gX-yCafi ~ 9axC^), (2.13) 
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r^ 7 = \g^(d a g^ + d igail -d^ ai ). (2.14) 



The quantities V are called the Christoffel symbols of the metric g. The connection coefficients 

reduce to them in the natural frame. They are zero for an orthonormal frame. 

lr The interpretation of this condition is that the second covariant derivatives of scalar functions commute, 
namely 

Vadgf- V d a f = O. 
In particular in the natural frame wJL, is symmetric in f3 and 7. 
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2.4 Curvature 



2.4.1 Definition 

The non-commutativity of covariant derivatives is a geometric property of the metric. It signals 
its curvature. The Riemann curvature tensor is defined as an exterior 2-form with value a linear 
map in the tangent plane to V by the following identity: 

(V A V M - V m VaK = Rx^pvP, (2.15) 

which gives by straightforward identification, 

Rx,,% = dxufo - dpufy + w%fJ^ - ufyfy - ufoC^. (2.16) 

2.4.2 Symmetries and antisymmetries 

The Riemann tensor is a symmetric double 2-form: it is antisymmetric in its first two indices, 
and in its last two indices written in covariant form. It is invariant by the interchange of these 
two pairs. 



2.4.3 Bianchi identities 

It holds that (vanishing of the covariant differential of the curvature 2-form) : 

V a -R/3 7 ,A^ + VpRjaM + VjR a p,\n = 0. (2-17) 

2.4.4 Ricci tensor, scalar curvature, Einstein tensor 

The Ricci tensor is defined by 

RaP = R\a, X /3- (2-18) 



The scalar curvature is, 



The Einstein tensor is, 



R = 9 aP R a p- (2.19) 



S a p = Rap - \g a pR- (2.20) 
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2.4.5 Conservation identity 

Contracting the Bianchi identities gives that, 

VaRfr^ - VpRyu + V 7 /2^ = 0, (2.21) 
and a further contraction gives the following identity satisfied by the Einstein tensor: 

V a S a ? = 0. (2.22) 
This identity implies that the sources must satisfy the so-called conservation laws, 

V a T a/3 = 0. (2.23) 

3 (n+l)-splitting adapted to space slices 
3.1 Definitions 

We consider a spacetime with manifold V = A4 x M and hyperbolic metric g such that the 
submanifolds A4 t = M. x {t} are spacelike. We take a frame with space axis e$ tangent to 
the space slice Ait and time axis eo orthogonal to it. Such a frame is particularly adapted to 
the solution of the Cauchy problem and will be called a Cauchy adapted frame. The dual 
coframe is: 

P = dx i + pdb, (3.1) 

ff is a time dependent vector tangent to M-t called the shift. The 1-form 6° does not contain 
dx l . We choose: 

8° = dt. (3.2) 

The pfaffian derivatives (action of the vector basis e a ) with respect to the adapted coframe are, 

8 = d t - (3 j dj, di = d/dx\ with d t = d/dt. (3.3) 

In the coframe 9 a the metric reads, 

ds 2 = -N 2 (8 ) 2 + gijPP. (3.4) 

The function N is called the lapse. We shall assume throughout TV > and the space metric g 
induced by g on M. t properly Riemannian. An overbar will denote a spatial tensor or operator, 
i.e., a t-dependent tensor or operator on M. Note that = gij and <f J = g l K 
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3.2 Structure coefficients 

The structure coefficients of a frame of a Cauchy adapted frame are found to be, 

C k = ~ C % = (3.5) 
and all other structure coefficients are zero. 

3.3 Splitting of the connection 

We denote by V covariant derivatives in the space metric g. Using the general formulas ( p,12| ), 
(|2l^ ), ( |2l4l) we find that, 



and 



that is, 



u% jk = r *jfc = r ifc; (3-6) 
w'oo = NgVdjN, lo° 0i = to° i0 = N^diN, cu° 00 = N^d N, 



u% = \N- 2 {d m + g hj C& + g ih C%), (3.7) 



<A = \N- 2 {d mj - g hj di(3 h - g ih d J f3 h ). (3.8) 

Using the expression (|3.3| ) of do we obtain that, 

= ^N~%g ih (3.9) 

where the operator do is defined on any f-dependent space tensor T by the formula, 

Bo = ^ - L p , (3.10) 

where Lp is the Lie derivative on A4t with respect to the spatial vector (3. Note that do T is 
a t-dependent space tensor of the same type as T. 

The extrinsic curvature Kij (second fundamental tensor) of Mt is classically defined as the 
projection on M. t of the covariant derivative of the unit normal v, past-oriented, that is 

K ij = ^{V l u J +V j u i ), (3.11) 

with, in our coframe, 

Vi = 0, uo = N. (3.12) 
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Lemma 3.1 The following identity holds, 

d g tj = -2NKij. (3.13) 

Proof. It holds that, 

V i u j = -uP ijVQ = --N^dogij = Kij. (3.14) 
The remaining connection coefficients are found to be, 

u'jo = —NlCj + djp, u\j = —NK l j, (3.15) 
and this completes the proof of the Lemma. ■ 

3.4 Splitting of the Riemann tensor 

We deduce from the general formula giving the Riemann tensor and the splitting of the con- 
nection the following identities, 



Rij,ki — Rij,ki + K ik Kij - KnK k j, 

Roijk = NiVjKki - VkKji), 

Roifij = N(d K i:j + NK lk K k j + VidjN). 
From these formulae one obtains the following ones for the Ricci curvature: 
NRij = NRij - doKij + NKijKl - 2NK tk K k j - V t d 3 N, 

N- 1 ^ = 3 d K h h - V h K h j, 

Roo = N(d K% - NKijRV + AN). 

Also, 

gVRij = R- N- 1 ^ + {Klf - N^AN, 



Soo = Roo — -^9ooR = 2 (-^oo + 9 %J 'Rij), 



hence, 



with R = gVRij 



2N~ 2 S 00 = -2S U = R- KijRV + (K%) 2 , 



3.16) 
3.17) 
3.18) 

3.19) 
3.20) 
3.21) 

3.22) 
3.23) 

3.24) 



4 Constraints and evolution 



We see in the above decomposition of the Ricci tensor that none of the Einstein equations 
contains the time derivatives of the lapse N and shift (5. One is thus led to consider the 
Einstein equations as a dynamical system for the two fundamental forms g and K of the space 
slices Ait- This dynamical system splits as follows. 
Constraints 

The restriction to Ait of the right hand side of the identities (3.19) and ( |3,23j ) contains only 



the metric gij and the extrinsic curvature Km of Ait as tensor fields on Ait- When the Einstein 
equations are satisfied, i.e., when, 

S aj 3 = R a p - -g al 3 R = T a p = p a/3 - -g a /3P, (4.1) 

these identities lead to the following equations called constraints: 
Momentum constraint 



Ci = ^(R 0i - Poi ) = -Vyrf + ViK h h - N- 1 ^, (4.2) 
Hamiltonian constraint 



C = ^ (Soo - Too) =R- K)K{ + {Klf + 2T °. (4.3) 

These constraints are transformed into a system of elliptic equations on each submanifold 
Ait, in particular on Aio, for g = g , K = Ko, by the conformal method (cf. C-B and York 
1980, Isenberg 1995, C-B, Isenberg and York 2000). 

Evolution 

The equations, 

Rij = R i:j — 2K jh Ki + KijK h — — - p {j , 

together with the definition, 

doQij = -IN Kij, 

determine the derivatives transversal to Ait of g and K when these tensors are known on Ait 
as well as the lapse N and shift 0, and source p^. It is natural to look at these equations as 
evolution equations determining g and K, while N and f3, projections of the tangent to the 
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time line respectively on eo and the tangent space to A4, are considered as gauge variables. 
This point of view is supported by the following theorem (Anderson and York 1997, previously 
given for sources in C-B and Noutchegueme 1988): 

Theorem 4.1 When Rij — p^ = the constraints satisfy a linear homogeneous first order 
symmetric hyperbolic system, they are satisfied if satisfied initially. 

Proof. When — p^ = we have, in the Cauchy adapted frame, with p = g a ^p a p^ 

R-p=-N 2 (R 00 -p 00 ), 

hence, 

s oo _ T oo = 1 (jR oo _ p oo } and R _ p = _ 2iV 2 (5 oo _ T o 0) = 2{s o _ T o ); 

and 

S» - T» = -\g l] {R -p) = -g t] (S° - r °), 

the Bianchi identities give therefore a linear homogeneous system for S = Sq — Tq and Eq = 
Sq — Tq with principal parts, 

N~ 2 d Zi + g^dj^, and SbSg + 

This system is symmetrizable hyperbolic, it has a unique solution, zero if the initial values are 
zero. The characteristic which determines the domain of dependence is the light cone. ■ 

5 Analytic Cauchy problem 

Geometrical initial data for the Einstein equations are a triple (A4,go,Ko) with go a properly 
Riemannian metric on the n-dimensional manifold A4, and Kq a symmetric 2-tensor on A4. 
A solution of the Cauchy problem for the Einstein equations with these initial data is an 
(n+l)-dimensional pseudo-Riemannian manifold (V, g) which we shall suppose of signature 
(—,+,...,+), such that M. can be identified with a submanifolds M.q of V, with go the metric 
induced by g on A^o and Kq the extrinsic curvature of Aio as submanifold of (V,g). The 
manifold (V, g) is called an Einsteinian development of the data. When the sources are given, 
for example zero (vacuum case), these geometrical initial data cannot be chosen arbitrarily, 
they must satisfy on M. the constraints (|4.2|), (|4.3|) . 
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The evolution equations read: 



dtgtj = ^NKij + ViPj + VjPi, 
d t k i:j = N{R lj -2K th K ] ^ + K lj K , f l }-V i d j N 
+ f3 h VK tJ + k ih V 3 f3 h + k hj ViP h - N Pij . 



(5.1) 



(5.2) 



No equation contains the time derivatives of the lapse iV and shift (3. We suppose that these 
quantities are given on V. The system is of the Cauchy-Kovalevski type, therefore we have the 
following theorem: 

Theorem 5.1 If the initial data are analytic on Mq while the sources, the shift and the lapse 
are analytic in a neighborhood of Ado then there exists a neighborhood of A4q inMxR such 
that the evolution equations have a solution in this neighborhood taking these Cauchy data. 

We deduce from this Theorem the following one: 

Theorem 5.2 If the sources satisfy the conservation laws, in particular are zero, on V the 
solution of the evolution equations satisfies the full Einstein equations if the initial data satisfy 
the constraints. 

Since the conservation laws depend also on the metric the application of these theorems 
requires further study, except in the case where the sources are zero. In this vacuum case we 
can enunciate: 

Corollary 5.3 Analytic initial data satisfying the vacuum constraints admit a vacuum Ein- 
steinian development. 

Proof. Take an arbitrary analytic shift (5 and N > 0. ■ 
We will come back later to the geometric uniqueness problem. 

6 Non-strict hyperbolicity of Rij = 

An evolution part of Einstein equations should exhibit causal propagation, i.e., with domain 
of dependence determined by the light cone of the spacetime metric. 

The equations Rij = are, when N and (5 are known, a second order differential system for 
gij. The hyperbolicity of a quasi-linear system is defined through the linear differential operator 
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obtained by replacing in the coefficients the unknown by given values. In our case and for given 
N,f3 and gij, the principal part of this operator acting on a symmetric 2-tensor 7^ is, 

\{{N- 2 d 2 m - g hk d 2 hk ) 7ij + d k d jlik + d k d l7jk - g hk d t d j7hk }. 

The characteristic matrix M at a point of spacetime is the linear operator obtained by replacing 
the derivation d by a covariant vector £. The characteristic determinant is the determinant of 
this linear operator. For simplicity we compute it in the classical case of space dimension n = 3. 
We find, 

DetM = $(g afi Uf>} 3 - 

The characteristic cone is the dual of the cone defined in the cotangent plane by annulation 
of the characteristic polynomial. For our system the characteristic cone splits into the light 
cone of the given spacetime metric and the normal to its space slice. Since these characteristics 
appear as multiple and the system is non-diagonal, it is not hyperbolic in the usual sense. One 
can prove by diagonalization of the system the following theorem (C-B 2000). 

Theorem 6.1 When N > and (3 are given, arbitrary, the system Rij = is a non-strict 
hyperbolic system in the sense of Leray-Ohya for g^, in the Gevrey class 7 = 2, as long as g^ 
is properly Riemannian. If the Cauchy data as well as N and (3 are in such a Gevrey class, 
the Cauchy problem has a local in time solution with domain of dependence determined by the 
light cone. 

7 Wave equation for K, Hyperbolic system 

Various hyperbolic systems have been obtained in recent years for the evolution of the dynam- 
ical variables (gij,Kij) by linear combination of with the constraints. The first of these 
hyperbolic systems has been obtained in C.B and Ruggeri 1983, zero shift, extended in C.B 
and York 1995 to an arbitrary shift. It uses a quasi-diagonal system of wave equations for Kij, 
shown to hold modulo a gauge condition for N, called now 'densitizing the lapse'. It works as 
follows. 

We use the expressions for Rqi and Rij, together with dogi k = —2NK ik which imply, 

d t% = -g hk {Vi{NK jk + Vj(NK ik ) - V fe (AT^)}, (7.1) 
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to obtain the identity, with fuj\ = + fji, and H = tv(K) = K%, 

Qij = doR^-V^o^-doiN-'doK^ + VhV^NK^-doiN^VjdiN) 

- NVid 3 H + BoiHKij - 2K im Kf) - V {i (K j)h d h N) - 2NR% m K^ 

- NR^Kfi + HVjdiN. (7.2) 

This identity shows that for a solution of the Einstein equations, 

the extrinsic curvature K satisfies a second order differential system which is quasi-diagonal 
with principal part the wave operator, except for the terms —NVidjH. The unknown g appears 
at second order, as well as N except for the term —do(N~ 1 VjdiN). It holds that, 

BoiN^VjdiN) + NVjdiH = N^VjdiidoN + N 2 H) + X ijt (7.3) 

where X{j is only of first order in g and K. 
We then densitize the lapse, that is we set, 

N = a(detg) 1/2 , (7.4) 

where a is an arbitrary positive tensor density of weight —1 called the 'densitized lapse'. Then, 



using the formula for the derivative of a determinant and the identity from Lemma 3.1, we find 
that: 

8 N = —N 2 H + 5 «(det g) 1/2 . (7.5) 

We see that d$N + N 2 H is an algebraic function of g. Hence, we have proved the following 
theorem. 

Theorem 7.1 Let a > and (3 be an arbitrary space tensor density and space tensor, depend- 
ing on t. Set N = a(det g) 1 ^ 2 , then: 

1. The equations 

fiii = 0, (7.6) 

are a quasi- diagonal system of wave equations for K. 

2. The equations above together with 

d gij = -2NK tj , (7.7) 
are a hyperbolic Leray system for g and K. 
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Proof. Part 1 has already been proved. To prove part 2, we give to the equations and 
unknown the following weights: 



m(K) = 2, m{g) = 2, ri(fi = 0) = 0, n(dg = 2NK) = 1. 



(7.8) 



The principal operator is then a matrix diagonal by blocks. Each block, corresponding to a 
pair (ij) of indices, is given by, 



with X a second order operator. The characteristic determinant is (— iV 2 ^ + <? u £i£j)£oi it 
is a hyperbolic polynomial. ■ 

8 Hyperbolic-elliptic system. 

An alternative method to the densitization of the lapse is to consider H as a given function h 
on space time, i.e., imposing given mean extrinsic curvature on the space slices. The second 
order equation for K obtained above reduces again to a quasi-diagonal system with principal 
part the wave operator. This gauge condition was used by Christodoulou and Klainerman, 
with h = 0, in the asymptotically Euclidean case, in the general case by C.B and York 1996. 
The lapse N is then determined through the equation Rq = p[] which now reads in the general 
(non-vacuum) case, 



This equation is an elliptic equation for N when g , K and p are known. 

Note that for energy sources satisfying the energy condition we have — p® > as well as 
| K | 2 = KijK l i > 0, an important property for the solution of the elliptic equation. The mixed 
hyperbolic-elliptic system that we have constructed will determine the unknowns N and g in 
a neighborhood of M. in M. x R when the shift (3 is chosen. 

9 Local existence and uniqueness 

Known properties of hyperbolic systems on manifolds and the use of the Bianchi identities to 
show the preservation of constraints under evolution, lead to the following theorem. 

Theorem 9.1 Let (A4,%,Kq) be an initial data set satisfying the vacuum constraints, where 
M. is a smooth n- dimensional manifold endowed with a smooth, Sobolev regular, Riemannian 




VdiN - (KijKV - p° )N = -doh. 
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metric e, and where go G , Kq G Hf_° c are a properly Riemannian metric and a symmetric 

2-tensor on A4 respectively. Suppose arbitrarily given on M. x X, X an interval o/R, the gauge 
variables a,0 G C°(I,Hf loc ) n C l {I,H^f) n C 2 {I,H^ C ). Then if s > § + 1 i/iere exists an 
interval J C T and a Lorentzian metric, 

g = -N 2 dt 2 + g lj {dx i + f3 i dt){dx j +f3 j dt), N = (a detg) 1/2 , (9.1) 

solution of the given Cauchy problem on A4 to , t G J , for the vacuum Einstein equations. For 
a given pair a, f3 this solution is unique. 

Notation 9.2 The spaces Hg ,loc are spaces of tensors which have generalized covariant deriva- 
tives in the metric e of order up to s which are square integrable on each open set of some given 
covering of A4, with norm uniformly bounded (i.e., for a given tensor, independent of the sub- 
set). The manifold (A4,e) is called Sobolev regular if the covering can be chosen so that these 
Sobolev spaces satisfy the usual embedding and multiplication properties. This is always the 
case if (A4,e) is complete. 

Two isometric space times (V,g) and (V,g) are considered as identical. One can prove the 
following theorem. 

Theorem 9.3 (Physical uniqueness) 1. Let (M.xJ'^g) be another solution of the Cauchy 
problem of Theorem \9. 1\ , with a different gauge choice. There exists an isometry of 
(M x J',g) onto (M x J',g), with J' C J and J' C J. 

2. The solution of the Cauchy problem for the vacuum Einstein equations is geometrically 
unique (i.e., up to isometrics) in the class of globally hyperbolic space times. 

For the definition of globally hyperbolic spacetimes see S. Cotsakis' lectures, this volume. 

10 First order hyperbolic systems 

Such systems are supposed to be more amenable to numerical computation. 
10.1 FOSH systems 

A first order system of N equations which reads, 

M aIJ (u)d aUl + f J {u) = 0, (10.1) 
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where Ui, I = 1, . . . , N, are a set of unknowns (for instance the components of one or several 
tensors), is called symmetric if the matrices M a are symmetric. Such a symmetric system is 
hyperbolic for the space slices Ait if the matrix M , coefficient of d/dt is positive definite. The 
energy associated to such systems is straightforward to write in a Cauchy adapted frame. In 



this case M l = M°, while M i = M i - ffM Q . Equations (|To7T|) imply, 



ujM aIJ (u)d aUl + ujf J (u) = ^d a (M aIJ u lUj ) + F(u) = 0. (10.2) 

We express the d^s as linear combinations of the usual time and space partial derivatives. 
Then, Eq. ( ^02| ) takes the form, 

ujM aIJ (u)— Ul + ujf J (u) = - — (M aIJ Ul uj) + F{u) = 0. (10.3) 

Integrating this equation on a strip A4 x [0, T] leads to an energy equality and an energy 
inequality (see later 'Bel-Robinson energy'). 

Anderson, C-B and York 1995 have written the system ( |7.6| ) , ( |7.7j ) together with the gauge 
condition (|7.5|) as a FOSH system, using also the equation Rqq = 0. 

A FOSH system had also been obtained by Frittelli and Reula 1994 just by combination of 
Rij with the constraints and densitization of the lapse. They used it to discuss the Newtonian 
approximation. A number of variants has been written since then, and their quality for numer- 
ical computation discussed (see in particular the Einstein-Christoffel system of Anderson and 
York 1999). 

10.2 Other first order hyperbolic systems 

Another criterion than symmetry has been recently used to test hyperbolicity of first order 
systems: It is the number of linearly independent eigenvectors associated to a multiple char- 
acteristic. If this number is equal to the multiplicity, and modulo conditions of uniformity, 
the system is hyperbolic. The verification of such a property involves heavy computations. 
Kidder, Scheel and Teukolsky 2001 write a whole family of such systems. They introduce as 
new unknowns the partial derivatives dhgij- They show that densitization of the lapse is a 
necessary condition for the hyperbolicity of the obtained systems. They evolve some of them 
numerically in the case of a one-black- hole spacetime and discuss their accuracy, i.e., how well 
the constraints are preserved. 
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11 Bianchi-Einstein equations 



The Riemann tensor is the geometric object which intrinsically defines gravitational effects. 
The following results bring nothing new for the local existence and unicity theorems but they 
are useful for the global-in-time studies. 

11.1 Wave equation for the Riemann tensor 

The Riemann tensor of a pseudo-Riemannian metric, R a /3 : \^, is antisymmetric in its pair of 
first indices as well as in its pair of last indices. We call it a symmetric double 2- form because 
it possesses the symmetry, 

The Riemann tensor satisfies the Bianchi identities 

V Q i?/3 7i AAt + V 7 il a) g i A A i + VpRj^Xfj, = 0. (H-2) 
These identities imply by contraction, 

Using the symmetry ( |1 1 . 1| ) gives the identities: 

V a R a PM + V„Rxp - VxR^p = 0. (11.4) 
If the Ricci tensor R a p satisfies the Einstein equations, 

then the previous identities imply the equations (Bel, Lichnerowicz) 

V a R a pM = VxPtf ~ y^Pxp- (H.6) 



Equations ([lL2|) and ( pT6|) are analogous to the Maxwell equations for the electromagnetic 
2-form F: 

dF = 0, SF = J, (11.7) 

where J is the electric current. 

Theorem 11.1 The Riemann tensor of an Einsteinian spacetime of arbitrary dimension sat- 
isfies a quasi-diagonal, semilinear system of wave equations. 
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Proof. One deduces from ( 11.2| ) and the Ricci identity, an identity of the form: 

where Spj^ is an homogeneous quadratic form in the Riemann tensor, 

)- (A -+ //)]}-{/? ^ 7 }- (H-9) 

Using equations ( |11.6| ), when the Ricci tensor satisfies the Einstein equations, gives equations 
of the form, 



with J/3 7) a^ depending on the sources p a p and being zero in vacuum: 

JfaM = V 7 (V M p A/ g - V A p M/3 ) - (J3 -> 7), 
and this completes the proof. ■ 



111.10) 



(11.11) 



11.2 Case n=3, FOS system 

In a coframe 9°,9 l where got = 0, equations ( |11.2 ) with {afi-y} = {ijk} and equations ( |11.6 ) 
with (3 = do not contain derivatives do of the Riemann tensor. We call them 'Bianchi 
constraints'. The remaining equations, called from here on 'Bianchi equations,' read as follows: 



Vn-R 



k-K0h,\fi, 



+ V h R, 



/i-ftfcO,A/i 



V i? i,xp + V h R iM = V A /y - Vppxi = J\ 



(11.12) 
(11.13) 



where the pair (\p) is either (Oj) or (J I), with j < I. There are 3 of one or the other of these 
pairs if the space dimension n is equal to 3. 

Equations (??) and (??) are, for each given pair (Oj), a first order system for the components 
Rhk,0j an d Rohfij- R w e choose at a point of the spacetime an orthonormal frame the principal 
operator is diagonal by blocks, each block corresponding to a choice of a pair (A/i, A < /x), is a 
symmetric 6 by 6 matrix which reads: 



/ 


do 








0-2 


-di 












do 








d 3 


-d 2 












do 


-d 3 





di 






0-2 





-d 3 


do 












-di 


d 3 








do 







V 





-d 2 


Oi 








do 


) 
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We have proved: 

Theorem 11.2 The Bianchi evolution equations are a FOS (first order symmetrizable) system. 
The Bianchi equations depend on the choice of frame, as does their hyperbolicity. 



11.3 Cauchy adapted frame 

The numerical valued matrix M f of coefficients of the operator d/dt corresponding to the 
Bianchi equations relative to the Cauchy adapted frame is proportional to the unit matrix, 
with coefficient N~ 2 , hence is positive definite and the following theorem holds. 

Theorem 11.3 The Bianchi equations associated to a Cauchy adapted frame are a FOSH 
system, with space sections Mf 

We will give an explicit expression of the full system after introducing two 'electric' and 
two 'magnetic' space tensors associated with the double 2-form R. They are the gravitational 
analogs of the electric and magnetic vectors associated with the electromagnetic 2-form F. That 
is, we define the 'electric' tensors by, 

E t , = R ifij , (11.14) 

Dij = \r, ihkr]jlm R hk ' lm , (11.15) 
while the 'magnetic' tensors are given by, 

Hij = \N-\ hk R hk fil , (11.16) 
Bji = ^N-\ hk R 0j hk . (11.17) 
In these formulae, n i j k is the volume form of the space metric g = gijdx % dxK 
Lemma 11.4 1. The electric and magnetic tensors are always such that 

Eij = Eji, Dij = Dji, Hij = Bji (11.18) 

2. If the Ricci tensor satisfies the vacuum Einstein equations with cosmological constant 

RaP = A<? a/3 (11.19) 

then the following additional properties hold 

Hij = Hji = Bij = Bji, E^ = (11.20) 



19 



Proof. (1) The Riemann tensor is a symmetric double 2-form, the electric and magnetic 
2-tensors associated to it by the relations possess obviously the given symmetries. 

(2) The Lanczos identity for a symmetric double two-form, with a tilde representing the 
spacetime double dual, gives 

Ra/3,\fi + Raf3,\fi = CaA^ _ Ca^9p\ + C/3/i3aA ~ C/3\9an, (11.21) 

with 

Caf3 = R a /3 — -^Rgaf3- (11.22) 

It implies that R a p,\^ + R a /3,x^ = if C a p = 0, in particular for an Einsteinian vacuum 
spacetime with possibly a cosmological constant. The relations ( |11.2C| ) can then be proved by 
a straightforward calculation that employs the relation ij 0i j k = Nr]^ between the spacetime 
and space volume forms. ■ 

In order to extend the treatment to the non-vacuum case and to avoid introducing un- 
physical characteristics in the solution of the Bianchi equations, we will keep as independent 
unknowns the four tensors E , D, B, and H, which will not be regarded necessarily as symmet- 
ric. The symmetries will be imposed eventually on the initial data and shown to be conserved 
by evolution. 

We now express the Bianchi equations in terms of the time-dependent space tensors E, H, 
D, and B. We use the following relations, found by inverting the definitions ( |11.14j )- (|11.17 ), 

Roifij = -N 2 Eij, Rhkfij = NifhkHij, (11.23) 

Rfik,lm = ff hkV 3 ImDij , Roj,lm = Nrfi m Bji. (11.24) 

We will express spacetime covariant derivatives of the Riemann tensor in terms of space covari- 
ant derivatives V and time derivatives, do, of E, H, D, B by using the connection coefficients 
in (3+l)-form as given in Section ||. 

The first Bianchi equation with [Xfj] = [Oj] has the form, 

Vo-R/ifc,Oj + Vfci?o/i,Oi - ^hRokfij = 0. (11.25) 

A calculation incorporating previous definitions, then grouping derivatives using do and Vj, 
gives to the first pair of Bianchi equations, with [Xfi] = [Oj], the following forms: 

BoEij - NrFiVhHy + (L 2 )y = J 0ji , (11.26) 
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where J is zero in vacuum and 



WhkHij) + 2NV [h E k]j + (LOftfcj = 0, (11.27) 

with, 

{L 2 )ij = -N{trK)Eij + NK k jE ik + 2NK k E k] (11.28) 

- (V h N)r, hl t H l3 + NK k hV lh iV m kj D lm + (V k N) V l kj B u , (11.29) 



{Li)hk,i = NK l 3 if hk H tl + 2(V [h N)E k]j + 2Nrf lj K\ k B h]i 



(11.30) 
(11.31) 



We see that the non-principal terms L\ and L 2 are linear in E, D, B, and H, with coefficients 
linear in the geometrical elements K and ViV. The characteristic matrix of the principal terms 
is symmetrizable. The unknowns E^ and with fixed j and i = 1,2,3 appear only in the 

equations with given j. The other unknowns appear in non-principal terms. The characteristic 
matrix is composed of three blocks around the diagonal, each corresponding to one given j. 

The j th block of the characteristic matrix in an orthonormal frame for the space metric g, 
with unknowns listed horizontally and equations listed vertically (j is suppressed) , is given by, 



(11.32) 



V 



e 











m 3 


-m 2 





e 


















e 


N£ 2 














Co 


















e 





-m 2 













Co 



This matrix is symmetric and its determinant is the characteristic polynomial of the E, H 
system. It is given by, 



a\2 



(11.33) 



The characteristic matrix is symmetric in an orthonormal space frame and the coefficient 
matrix M f is positive definite (it is the unit matrix). Therefore, the first order system is 
symmetrizable hyperbolic with respect to the space sections Mt- We do not have to compute 
the symmetrized form explicitly because one can obtain energy estimates directly by using the 
contravariant associates E lJ ' , H l i , ... of the unknowns. 
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The second pair of Bianchi equations, for D^j and B^j, obtained for [Xfi] = [Im] is analogous. 
The characteristic matrix for the [Im] equations, with unknowns Dij and Bij, j fixed, with an 
orthonormal space frame, is the same as the matrix found above. 

If the spacetime metric g is considered as given, as well as the sources, the Bianchi equations 
form a linear symmetric hyperbolic system with domain of dependence determined by the 
light cone of g. The coefficients of the terms of order zero are ViV or iVK. The system is 
homogeneous in vacuum (zero sources). 

11.4 FOSH system for u = (E, H, D, B, g, K, f ) 

The Bianchi equations depend on the metric. Our problem is to find a system for determining 
the metric from the Riemann tensor (through eventually other auxiliary unknowns), which 
together with the Bianchi equations, constitute a well-posed system. It is possible to construct 
a FOSH system linking the metric and the connection to our Bianchi field (Anderson, C.B and 
York 1997), if we again densitize the lapse, i.e., set N = a(detg) 1//2 . This system is inspired 
by an analogous one constructed in conjunction with the Weyl tensor by H. Friedrich 1996. 

11.5 Elliptic - hyperbolic system 

Instead of determining the metric from the curvature through hyperbolic equations, one can 
try to do so by elliptic equations on space slices. Well-posed problems for such equations are 
essentially global ones and depend on the global geometric properties of the space manifolds. 

11.5.1 Determination of K 

We deduce from the identities ( |3.17| ), which read, 

N^Roijk = VjK ki - VkKji, (11.34) 

and the Ricci identity, that 

V^JV- 1 ^-*) = V j VjK ki - V^Kji - R j kjh K hl - R j kih Ky t . (11.35) 

As a gauge condition we now suppose that H = K% is a given function h on the spacetime, N 
then satisfies on each slice the equation deduced from the Einstein equation Rqq = 0: 

AN - K ij K i:j N = -d h (11.36) 
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The use of the momentum constraint on each slice gives for K (when N and the Riemann 
tensor of spacetime are known) a quasi-diagonal semilinear system, elliptic if g is properly 
Riemannian, namely, 

V J V 3 K ki - W m K hi - SP^Kji = VAh + V 3 {N- l R 0i)jk ), (11.37) 

where (cf. flOp), 

RijM = Rij,kl — K ik Kij + K a K k j. (11.38) 

The global solvability of the equation for K can be proved under some conditions, for instance 
in a neighborhood of Euclidean space, perhaps also for a Robertson- Walker spacetime with 
compact space of negative curvature. 

11.5.2 Determination of g 

The equation dog™ = —2NKij determines as before g when N,(3 and K are known. However 
it does not improve the regularity on M. t of g over the regularity of K. 

A better result can be sought through the identity which gives the Ricci tensor of g in terms 
of the Riemann tensor of spacetime and K by, 

Rij = 9 hk Rih,jk + HK l3 - K^Kjh- (11.39) 

Various methods have been devised to determine a Riemannian metric from its Ricci tensor by 
elliptic equations (see in particular Andersson and Moncrief, to appear, for the case of compact 
manifolds with negative curvature). 

12 Bel- Robinson energy 

The new formulation brings nothing new for the local existence and uniqueness theorem. It 
is useful in obtaining geometrical energy estimates (Bel-Robinson energy) leading possibly to 
global existence theorems. Such estimates have been used by Christodoulou and Klainerman 
1989 in the case of asymptotically Euclidean manifolds. They are used by Andersson and 
Moncrief (to appear) for compact manifolds with negative curvature. 
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12.1 Bel-Robinson energy in a strip 

Multiply ( FX27p by \m hk H l i, and recall that rjfrfhk = 2(5*/ , r/^r/^ = <$V fc r - 5 i r 5 h h and 



005^ = 2NK l K Then we find that, 



^HVdoWhkHy) = ^BoiHijHV) - M U (12.1) 



M x = jr j l rs H lm r, i hk H ij d (g hr g ks g jm ) (12.2) 



= N(K%H ij - K\H lj + Ki> H il )Hij. 
Likewise, multiply ( 11.26| ) by E l i to obtain, 



E^doEij = -daiEijE 1 ') - M 2 , (12.3) 

M 2 = N(K\E lj + K l j E il )E ij . (12.4) 

Multiplication by appropriate factors (cf. Anderson, C.B and York 1997) of the second pair of 
Bianchi equations leads to analogous results. The sum of the expressions so obtained from the 
four Bianchi equations gives an expression where the spatial derivatives add to form an exact 
spatial divergence, just as for all symmetric systems. Indeed, we obtain, 

\ d (|E| 2 + |H| 2 + |D| 2 + |B| 2 ) + VhiNEV^iHu) 

- V h (NB i i V lh i D lj ) = Q(E,H,T>,B) + S, (12.5) 

where we have denoted by | • | the pointwise g norm of a space tensor, and where Q is a quadratic 
form with coefficients ViV and K. The source term iS, zero in vacuum, is, 

S = JoijEP - ^NJ lmiVh lm B ih . (12.6) 

We define the Bel-Robinson energy at time t of the field (E, H, D,B), called a 'Bianchi field' 
when it satisfies the Bianchi equations, to be the integral, 



B(t) = - I (|E| 2 + |H| 2 + |D| 2 + |B| 2 W t . (12.7) 



gf 
M t 



We will prove the following. 



Theorem 12.1 Suppose that g is C 1 on M x [0, T\ and that the gt norms of ViV and K are 

uniformly bounded on Mt, t £ [0, T]. Denote by ir(t) the supremum 

7T(t) = Sn VMt (\VN\ + \K\). (12.8) 
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Suppose the matter source J G L l ([0, T], L 2 (ftA t )), then the Bel energy of a C 1 Bianchi field 
with compact support in space satisfies for < t < T the following inequality, 

Bit) 1 / 2 < (B(0) l / 2 + % f\\J\\ L 2 {MT) dT)zMC fv{T)dT), (12.9) 
1 Jo Jo 

where C is a given positive number. 



Proof. We integrate the identity ( 12.5| ) above on the strip M. x [0, t] with respect to the 
volume element [ig T dr. If the Bianchi field has support compact in space the integral of the 
space divergence term vanishes. The integration of dof on a strip of spacetime with respect to 
the volume form dr \i^ t goes as follows, for an arbitrary function /, 

/ d fn gT dT=[ [ (d t - pdi)f n gT dr. 

J Mr " Jo J Mr 



It holds that, 



dof /i 



/ {dt(fv St : 

JMt 



gf 



(12.10) 



>Mt JMt 

Using the expression for the derivative of a determinant and the relation between g and K, we 
find that, 



d tl i- gt = (-JVtrK + V^)^- gt 
and therefore if / has compact support in space, 



/ d f^ t = d t [ f» St + f /iVtrK Mgt . 

JMt JMt JMt 



The integration on a strip leads therefore to the equality, 



B(t) = B(0) + 



JMt 



(Q + S) Hg T dr, 



(12.11) 



(12.12) 



(12.13) 



withg, 



Q = Q + -iV(trK)(|E| 2 + |H| 2 + |D| 2 + |B| 2 ) 



We deduce from this equality, and the expression indicated in (12.5) for Q, the following 
inequality, with C some number, 

B(t) < B(0) + C{ /"sup MT (|ViV| + |K|)B(r) + || J\\ mMT) B{t) 1 / 2 } dr. (12.14) 
Jo 

This inequality and the resolution of the corresponding equality imply the result. ■ 

2 One can take advantage of the decomposition of Q to obtain better estimates in the case TrK< 0. 
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Remark 12.2 The quantities — V^iV and 2K^ are respectively the (Ok) and (ij) components 
of the Lie derivative C n g of the spacetirae metric g with respect to the unit normal n to Mt 
(its (00) component is identically zero). The Bel-Robinson energy is therefore conserved if this 
Lie derivative is zero. 

The estimate of the Bel-Robinson energy is only an intermediate step in global existence 
proofs since it depends on the metric which itself depends on curvature. 

12.2 Local energy estimate 

We take as a domain Q of spacetime the closure of a connected open set whose boundary <9f2 
consists of three parts: A domain uo t of Mt, a domain loq of Mo, and a lateral boundary C. 
We assume C is spacelike or null and 'ingoing', that is timelike lines entering 17 at a point of 
C are past-directed. We also assume that the boundary dO, is regular in the sense of Stokes 
formula. We use the identity previously found and integrate this identity on Q, with respect to 
the volume form \x^ t dr. It can be proved that the integral on C resulting from the application 
of Stokes formula is nonnegative. The Bel-Robinson energy on oo t satisfies therefore the same 
type of inequality as found before on Mt- In particular, we have B(u T ) = 0, if B(ujq) = and 
J = (vacuum case). Then E = H = D = B = 0inr2if they vanish on the intersection of 
Mq with the past of Q (result found by York 1987). Note that such a result is not sufficient 
to prove the propagation of gravitation with the speed of light because it treats only curvature 
tensors that are zero in some domain, not the difference of nonzero curvature tensors. The 
Bianchi equations are not by themselves sufficient to estimate such differences because their 
coefficients depend on the metric, which itself depends on the curvature. 



13 (n+l)-splitting in a time-adapted frame 
13.1 Metric and coframe 

We choose the time axis to be tangent to the time lines, i.e., the cobasis 6 is such that 9 l does 
not contain dx°. We set, 



9 l = a)dx 3 , 



f = Udx + bidx\ 



(13.1) 



We will call such a frame a CF- (Cattaneo-Ferrarese) frame. The Pfaff derivatives d a in the 
CF-frame are linked to the partial derivatives d/dx a by the relations, 
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with A] the matrix inverse of a] . The structure coefficients of the coframe are found to be, 

c 0i = N- 1 d i N-Ajd b j = -c%, 

and with f {ij] = fa - f jh 

4 = UA^iN^h), 4 = Aidoa), c\ k = A\ h d k] a). 

Remark 13.1 If the time lines are not hypersurface orthogonal (i.e., if hi ^ 0), the coefficients 
c\j are different from the structure coefficients of the space frame 9 l . 

Choosing the frame to be orthonormal the metric reads, 

3 

g = -(8P) 2 + Y,(8 i ) 2 - (13-2) 



i=i 



13.2 Splitting of connection 

We deduce from the general formulas, 



, ,o _ , ,o _ n 



Yi = uj 00 ,i = = ujq = -coifl = c° 0i = U 1 d i U - Ajdobj, 
and we know that lu^ = uj i,j is antisymmetric in i and j. We set, 

wwj = fij = \{Apoal - A%a> + A^b h }. 



Let e a = d a be the frame dual to 9 a , i.e., such that the vector e a has components <5^ Q ^. Then, 

V /3 e a = Ufa, 

in particular, ujqi j is the projection on erj\ of the derivative of eu\ in the direction of . We 
have fij = if the frame is Fermi-transported along the time line. We will make this hypothesis 
to simplify the formulas. 

The connection coefficient uJio,j is the sum of a term symmetric in i and j and an antisym- 
metric one and therefore we have, 

Xij = u i0 j = coj = ^{Apoai + A%4 + A^b h }. 

The antisymmetric term vanishes if the time lines are hypersurface orthogonal (pi = 0). 

The coefficients u>^ are also linear expressions in terms of the first derivatives of the frame 
coefficients, they are identical to the connection constructed with the a\ at fixed t only if bi = 0. 
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13.3 Splitting of curvature 

Using the general formulas we find in the chosen frame, 

Rol.j = <>oAj + XnMj + Y'Xhj ~ YjXh' (13.3) 

we denote by V is the pseudo-covariant derivative constructed with di and (Cataneo- 
Ferrarese transversal derivative). We have, 

/>%,:,, = V,.V ; - <V .V r ; - X^Xf (13.4) 

R hk...j = R'hk...j + x '.'.'.l x jh - Xj k X l h , (13.5) 

where R^ k ^ denotes the expression formally constructed as a Riemann tensor with the coef- 
ficients and, 

R'kh...j = ^kXhj - VhXkj - Yj(X kh - X hk ). (13.6) 

Remark 13.2 The symmetry Rkhflj = Roj,kh results from the expression of the connection in 
terms of frame coefficients. 

We deduce from the splitting of the Riemann tensor the following identities: 

R 00 = ViY 1 - d Xf - X h j Xj- \ (13.7) 

R h0 = VjX^ j - V;,.v; - Y\X 3h - X hj ). 

13.4 Bianchi equations (case n=3) 
13.4.1 Bianchi quasi constraints 

The Bianchi identities and their contraction contain, as in a Cauchy adapted frame, equations 
which do not contain the derivative do of the Riemann tensor, namely, 

ViRjhM + V h Ri jM + VjR hiM = 0, (13.8) 

V a R\x» = -V^ A0 + V A/> . (13.9) 
We call these equations Bianchi quasi-constraints. 
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13.4.2 Bianchi evolution system 

The remaining Bianchi equations can be written, as in the case of a Cauchy adapted frame, as 
a FOS (first order symmetric) system for two pairs of 'electric' and 'magnetic' 2-tensors. This 
system cannot be said to be hyperbolic in the usual sense: The principal matrix M° is the unit 
matrix hence positive definite but the operators d/dt appears also in the matrices M*. We say 
that the system is a quasi-FOSH system. It is a usual FOSH system with t as a time variable 
if the matrix of coefficients of d/dt is positive definite. It can be proved that this is the case if 
the metric induced on the t = constant submanifolds, 



is positive definite and U > 0. 

13.4.3 Quasi-FOSH system for connection and frame 

When the Riemann tensor is known the identities which express it become equations for the 
connection. Some of them do not contain the derivative do, we call them connection quasi- 
constraints. Identities linking connection and frame become first order equations for the frame 
coefficients. No equation gives the evolution of U. It can be considered as a gauge variable 
fixing the time parameter. 

13.5 Vacuum case 

In vacuum we give arbitrarily on the spacetime V the scalar U, length of the tangent vector d/di 
to the time line, together with the projection Y« of V eo eo on ej. The quantities // being chosen 
zero, the identities previously written give, when the Riemann tensor is known, equations with 
principal operator the dragging along the time lines of uj l h - and Xy, and when the connection is 
known equations for the dragging of the frame coefficients. These equations together with the 
Bianchi evolution equations constitute a quasi-FOSH system. This system is a FOSH system 
with respect to t as long as gij is positive definite. 

13.6 Perfect fluid 

In the presence of fluid sources one can obtain a quasi-FOSH system for the gravitational and 
fluid variables by taking as time lines the flow lines and proceeding as follows (Friedrich 1998). 




(13.10) 



h 



29 



13.6.1 Fluid equations 

The stress energy tensor of a perfect fluid is, 

T a/ 3 = + p)u a up + pg a /3. 

Then, 

and one supposes that the matter energy density \i is a given function of the pressure p. 

The Euler equations of the fluid, which express the generalized conservation law V a T alS = 0, 
are equivalent to the equations, 

(p + p^V^ + (^u 13 + g a/3 )d a p = 0, with u a u a = -1, 

and, 

(fi + p)V a u a + u a d a fi = 0. 

In our coframe they read, 

(ji + p)Yi + dip = 0, Yi = ui , 

don + (M + P) x t = 0- (13.11) 
Using the index F of the fluid defined by, 

dp 



F(p) = J 



nip) +p 

we have, 

Yt = -diF, (13.12) 

and, 

UpdoF + Xj = 0. (13.13) 
The commutation relation between Pfaff derivatives and the definitions give that, 

id di - did )F = d&dcF = YdoF - XjdjF, 



30 



and therefore, 

fi' p [d Yi + Yid F + (fl - xi)djF] - din' p d F - d t Xl = 0. (13.14) 

The use of previous identities replaces d a F by functions of Y, X and p. The derivatives difi' p 
are functions of Y and p since. 

Following H. Friedrich, we replace diX^ by its expression deduced from the equation, 

Rio = ^hXl" h — VjX^ - - Y h (X hi - X ih ), 
and changing names of indices we obtain, 

fi' p d Y h - VjX J h - Y h {X hi - X ih ) + Y h d F + 

-XldjF) + d hf i' p d F = 0. (13.15) 
Replacing V h Yi by ViY h + c° hi Y where, 

Y = -d F = -(/x;)- x x|, 

c ih = ^hi ~ ^ih = —^hifl + Uih,0 = X ih — X hi , (13.16) 

we obtain, 

doXf* - vV h + {» p )- x Xt{X ih - X hi ) - Y h Y i + X-^Xf + X- h - j = A';;,,;.. „. (13.17) 

The principal operator on the unknowns Y and X in the above equations is diagonal by blocks 
and symmetric. The /i-block reads: 

/ fi' p d -di -d 2 -d 3 \ 

-di d 

-d 2 d 

\ do J 

If /Lip > the matrix M° is positive definite in the CF-frame. The system is a quasi-FOSH 

j 
h- 



system for the pairs Yh , Xl . 



Remark 13.3 The characteristic determinant associated with the system ( 13.lt ), ( 13.17 ) is, 



i=l,2,3 



The roots of // p £ — Yli=i 2 3 £? = ® correspond to sound waves. Their speed is at most 1 (speed 
of light) if and only if fj,' > 1. 
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13.6.2 Sources of the Bianchi equations 

In our frame the source tensor p reduces to, 

Poo = + Poi = 0, Pij = \^~ P) S H- 

We have seen that d a p and d a fi are smooth functions of p, Y and X. The same property holds 
for Jifij and J^k- 

13.7 Conclusion 

Assembling the results of the previous subsections we find the following theorem. 

Theorem 13.4 The Einstein equations with source a perfect fluid give a quasi-FOSH system 
for the Riemann curvature tensor, the frame and connection coefficients, and the density of 
matter, when the flow lines are taken as timelines, U and /j = u> l j given arbitrarily. 

Corollary 13.5 The EEF ( Einstein- Euler-Friedrich) system is a FOSH system relatively to 
t = constant slices as long as the quadratic form, 

9jh= Yl a K _6 A, (13.18) 

i=l,2,3 

is positive definite, U > and fj,' p > 1. 
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